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We evaluate interhadron interactions in two-color lattice QCD from Bethe-Salpeter amplitudes 
on the Euclidean lattice. The simulations are performed in quenched SU(2) QCD with the plaque- 
tte gauge action at f3 = 2.45 and the Wilson quark action. We concentrate on S-wave scattering 
states of two scalar diquarks. Evaluating different flavor combinations with various quark masses, 
we try to find out the ingredients in hadronic interactions. Between two scalar diquarks (uC'ysd, the 
lightest baryon in SU(2) system), we observe repulsion in short-range region, even though present 
quark masses are not very light. We define and evaluate the "quark-exchange part" in the inter- 
action, which is induced by adding quark-exchange diagrams, or equivalently, by introducing Pauli 
blocking among some of quarks. The repulsive force in short-distance region arises only from the 
"quark-exchange part", and disappears when quark-exchange diagrams are omitted. We find that 
the strength of repulsion grows in light quark-mass regime and its quark-mass dependence is similar 
to or slightly stronger than that of the color-magnetic interaction by one-gluon-exchange (OGE) pro- 
cesses. It is qualitatively consistent with the constituent-quark model picture that a color-magnetic 
interaction among quarks is the origin of repulsion. We also find a universal long-range attractive 
force, which enters in any flavor channels of two scalar diquarks and whose interaction range and 
strength are quark-mass independent. The weak quark-mass dependence of interaction ranges in 
each component implies that meson- exchange contributions are small and subdominant, and the 
other contributions, ex. flavor exchange processes, color-Coulomb or color-magnetic interactions, 
are considered to be predominant, in the quark-mass range we evaluated. 

PACS numbers: 12.38. Gc, 12.39.Pn 



I. INTRODUCTION 

Hadron-hadron interactions play important roles in nu- 
clear and hadron physics. Among them baryon-baryon 
interaction is one of the most important issues to be clar- 
ified in reliable manners: Baryons are the main building 
blocks of our world and their interactions are essential 
for the structures and dynamics of nuclei. After the suc- 
cessful proposal by Yukawa in 1935 pj, the long- and the 
intermediate-range interactions have been described by 
one-boson-exchange potentials (OBEP) 0-01 • The key 
players for the hadronic interactions in these regions are 
considered as mesonic modes. On the other hand, the 
short-range part of the interactions has not been well 
clarified so far. Especially, the origin of the large repul- 
sive core in nucleon-nucleon (NN) channel Q, which ac- 
counts for the high-density nuclear phenomena, has been 
one of the longstanding problems in hadron physics 0,0- 
[l5| . The degrees of freedom of quarks and gluons would 
be essential in such short-distance region, and hence they 
should be clarified directly in terms of the fundamental 
theory, Quantum ChromoDynamics (QCD). 

It is now common knowledge that all the hadronic sys- 
tems are governed by QCD, and it is long after QCD 
was established as the fundamental theory. Nevertheless, 
QCD-based description of hadronic systems has not yet 
been successful due to the strong coupling nature of the 
low-energy QCD, whose dynamics is far from the pertur- 
bative regime and essentially nonperturbative. Among 
potentially workable strategies for nonperturbative anal- 



yses of QCD, lattice QCD is a unique reliable method 
for low-energy QCD. It can now reproduce emp irical 
hadronic masses with a very good accuracy [16[, and 
its outcome can be directly compared with experiments. 
There have been several attempts which aim at clari- 
fication of interhadron interactions by means of lattice 
QCD [i~7l - [2i"l ] . Recently, a nucleon-nucleon (NN) interac- 
tion was evaluated from NN Bethe-Salpeter (BS) ampli- 
tudes on the lattice [22| , where the short-range repulsion 
and the attraction in the intermediate region were ob- 
served. While the information of scattering phase shifts 
are properly encoded in asymptotic BS amplitudes, "po- 
tentials" directly constructed from BS amplitudes have 
setup- (operator- or energy-) dependences. Nevertheless, 
they are considered useful to gain qualitative understand- 
ing of hadronic interactions. 

In this article, we employ SU(2) lattice QCD, and in- 
vestigate hadronic interactions aiming at clarifying the 
essential structure of them. First, the scalar diquark 
(uC^d, the lightest baryon in SU(2) system) is an 
isospin-zero scalar hadronic state in SU(2) QCD and 
some of possible one-meson-exchange channels in diquark 
interactions are restricted, and hence analyses could be 
simpler. Diquarks in SU(2) QCD have a direct connec- 
tion to mesons, and two-diquark correlators are identical 
with two-meson correlators if we neglect possible quark- 
disconnected diagrams. Second, while SU(2) and SU(3) 
QCDs have different natures in some situations [HI, Hj], 
they often show similar aspects [25|, and hence SU(2) 
QCD can be a testbed for understanding hadron interac- 
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tions. On the lattice, it is possible to switch on/off Pauli- 
blocking effects among quarks controlling flavor contents 
or equivalcntly quark-exchange diagrams, which enables 
us to classify hadronic interactions in terms of their ori- 
gins. 

The paper is organized as follows. In Sec. [TTl we give 
a brief explanation of our strategy. Lattice QCD results 
are shown in Sec lIIIl and SeclIVlis devoted to discussions. 
In Sec|Vl we show the results of additional trial analyses. 
The summary and conclusions are given in Sec. IVII 



II. FORMULATION 

We follow the strategy proposed by CP-PACS 
group [26j . where asymptotic Bethe-Salpeter wavefunc- 
tions on the Euclidean lattice are adopted to evaluate 
pion scattering length. This method has been further 
employed for various hadron-hadron channels [27L [28| . 

We measure a two-diquark correlator Wjj^^R, T) at 
Euclidean-time t, which is defined as 

Wi JlW> r(R,t)= (1) 
E(%r(x, i)D fci , r (x + R, t)D^ r (0, 0)^(0, 0)02) 

X 

Here, R represents the relative coordinate of two scatter- 
ing hadrons, and (x, t) and £)J- r (x, t) are interpo- 
lating fields for diquark states whose flavors are i and j. 
Dij^rfct) is defined with two quark fields, 

A,-,r(x, t) = s ab qf(x, t)r<#(x, t), (3) 

with e ab being the 2x2 anti-symmetric tensor. Spinor 
indices of the two quark fields in D are contracted with 
T. Possible T's are C, C75, C-f^ and 0^75, which re- 
spectively correspond to pseudo-scalar, scalar, axialvec- 
tor, and vector diquarks. We employ wall-type opera- 
tors for sources, while we use point-type operators for 
sinks. After Euclidean-time evolution, a system is fi- 
nally dominated by the ground-state, and W(R, t) loses 
t-dependence besides an overall constant and an expo- 
nential damping factor. The BS wavefunction W(R) is 
then obtained as W(R) = lim^oo W(R,t). 

With such W(R), we can define an R-dcpcndcnt func- 
tion V(R) assuming a nonrclativistic Schrodinger-type 
equation, 

(C + ^ (R) ~ E ) W(R) = °' (4) 

Thus extracted scattering-energy dependent function 
V(R) is an equivalent potential which reproduces the 
scattering phase shift at the target energy. Though the 
information of phase shift is properly encoded in such 
equivalent potentials V(R), they could be energy- and 
operator-dependent. In this sense, equivalent potentials 
V(R) contain no more definite information than phase 



shifts in asymptotic region. Nevertheless, equivalent po- 
tentials are expected to be useful to gain understanding 
of hadronic interactions, at least qualitatively. We ex- 
tract equivalent potentials, which we simply call as "po- 
tential" in this paper, and evaluate hadronic interactions 
in SU(2) QCD. 

We concentrate on S-wave scattering states of two 
scalar diquarks (r = C75), and we project wavefunctions 
W(R) onto A£ - wavefunction W(R) = W(|R|), which 
has overlap with I = states, by summing over R in 
terms of corresponding discrete rotations. We here ne- 
glect the contributions from I > 4 scattering states, since 
such contributions can be dropped by taking large Eu- 
clidean time t. 

All the simulations are performed in SU(2) quenched 
QCD with the standard plaquette gauge action and the 
Wilson quark action. The lattice size is 24 3 x 64 at 
f3 = 2.45, whose lattice spacing is about 0.1 fm if we as- 
sume y/cf is 440 MeV [23, H(| ■ We employ four different 
Hopping parameters k = 0.1350, 0.1400, 0.1450, 0.1500 
for quarks. Diquarks at these k's would be well described 
by conventional quark models, and Nambu-Goldstonc- 
boson nature does not emerge. (See Sees . IIII A"l and IVBl ) 



III. LATTICE QCD RESULTS 

A. Hadron masses 

We show the lowest-state diquark masses in each 
channel in Table |TJ The masses are extracted in a 
standard manner. We fit hadronic correlators C(t) = 
(H(t)H'(t SIC )) by a single-exponential function, C(t) = 
C exp(— mt). We do not need "quark-annihilation" 
diagrams in the computation of diquark correlators. 
The scalar, vector, pseudoscalar, axialvector (r = 
C75j ^7^75, C, C7 M ) channels are investigated. 

Fig. [1] shows the mass difference Am between scalar 
and axialvector diquarks. The dashed line denotes the 
fit function, Am = Ctoq 2 , where iiiq = |mpv is 
the half of the axialvector-diquark (vector-meson) mass. 
The mass splitting Am is clearly proportional to niq 2 , 
which supports the color-magnetic interaction as the ori- 
gin of the mass splitting and indicates the validity of 
nonrelativistic-quark-model description of scalar and ax- 
ialvector diquarks, at least in the quark-mass range we 
consider. In the picture of constituent quark model, if 
spatial wave functions of constituent quarks in diquarks 
are affected by quark-mass variation, the Tiq 2 depen- 
dence of the mass splitting can be modified from the 
m,Q 2 scaling, which originates from the factor mg 2 of the 
strength of the color-magnetic interaction. The almost 
perfect mg 2 -behavior of the mass splitting then implies 
that the spatial wave functions in these two diquarks are 
similar and less dependent on quark mass. 
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FIG. 1: Scalar-axialvector mass splitting is plotted as a func- 
tion of axialvector-diquark mass. The dotted line is a fit func- 
tion CmT? . 



k Scalar Axialvector Pseudoscalar Vector 


Am 


0.1350 1.044(2) 1.056(2) 
0.1400 0.836(2) 0.855(2) 
0.1450 0.618(2) 0.651(2) 
0.1500 0.377(3) 0.447(2) 


1.285( 2) 
1.102( 2) 
0.919( 4) 
0.757( 7) 


1.286( 3) 
1.102( 4) 
0.918( 5) 
0.728( 4) 


0.012(2) 
0.019(2) 
0.033(2) 
0.070(3) 


0.1150 1.020(2) 1.026(2) 
0.1250 0.494(2) 0.504(1) 


1.397( 9) 
1.053(25) 


1.352( 8) 0.006(2) 
0.889(19) 0.010(2) 



TABLE I: All the hadronic masses are listed. The masses at 
k=0.1150 and 0.1250 are obtained without high-energy glu- 
ons. (See Sec. IV Bl ) Am represents the scalar-axialvector 
diquark mass splitting. 
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FIG. 2: The wavefunctions W(R,t) at k=0.1500 are plotted 
as a function of a source-sink separation t. 



B. Wavefunctions 

We hereby consider two-scalar-diquark wavefunctions 
with two different flavor combinations. One combination 
is (i,j,k,l) = (1,2,1,2), where only two independent 
flavors exist in four quarks. The other is (i,j,k,l) = 
(1,2,3,4), where all the quarks have different flavors. 
We adopt the same hopping parameters (quark masses) 
for all the quarks in both cases, and therefore all the 
scalar diquarks degenerate in mass. We note that the 
diagrams for the (i,j,k,l) = (1, 2, 1, 2) channel consist of 
direct and quark-exchange diagrams, and the direct dia- 
grams in the k, I) = (1, 2, 1, 2) channel are identical 
with those needed for the (i,j,k,l) = (1,2,3,4) chan- 
nel. "Quark-exchange diagrams" in quark-propagator 
contraction are included only in the k, I) = (1, 2, 1, 2) 
case. The existence of quark-exchange diagrams would be 
essential for the short-range interactions, since any Pauli- 
blocking effects among quarks are not included without 
exchange diagrams. We expect that the origin of short- 
range hadronic interactions can be accessed by comparing 
these two flavor combinations. 

We refer to the determination of wavefunctions W(R). 
Law wavefunctions W(R,t) depend on the source-sink 
separation t due to possible excited-state contamina- 
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FIG. 3: The potentials V(R,t) at «;=0.1500 and R=l, 2, 4 
are plotted as a function of a source-sink separation t. 



tions. We then have to take enough large t to ensure 
the absence of such contaminations. For this aim, we ex- 
tract t-dependent wavefunctions W(R, t) and determine 
t-window where excited-state contaminations are negli- 
gible. Fig. [2] shows t-dependent wavefunctions W(R, t) 
at k=0.1500 as functions of source-sink separation t. 
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The upper and lower panels correspond to (i,j,k,l) = 
(1,2,3,4) and (1,2,1,2) cases, respectively. They show 
plateaus at t > 25, and we then extract "wavefunctions" 
W(R) by fitting the data as W(R,t) = W{R) at t > 25. 
For further confirmation, we look at ^-dependent poten- 
tials V(R, t) extracted directly from W(R, t). In Fig. EH t- 
dependent potentials V(R, t) (i?=l,2,4) at k=0.1500 arc 
shown. One can find that they consistently show plateaus 
at t > 25. 
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hadron size itself and would indirectly affect interhadron 
potentials, these effects seem small in the present quark- 
mass range. 
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FIG. 5: Potentials Vi2,34(-R) computed with the flavor com- 
bination, (i,j,k,l) — (1,2,3,4), are plotted as functions of 
relative distance R. 



FIG. 4: The wavefunctions W(R) for the flavor combina- 
tion k, I) = (1, 2, l, 2) are plotted as functions of relative 
distance R for ^=0.1350, 0.1400, 0.1450 and 0.1500. The cor- 
responding hadronic masses can be found in Table [I] 

In Fig. EH as an example, we show the wavefunctions 
W(R,25) as functions of the relative distance R for four 
different K's.The flavor combinations are all (i,j,k,l) = 
(1,2,1,2). 



C. Potentials : (i,j,k,l) — (1,2,3,4) case 

We next proceed with the potentials extracted from the 
wavefunctions. Fig. [5] shows the reconstructed potentials 
plotted as functions of R, where the flavor combination 
is set to (i,j,k,l) = (1,2,3,4). In this case, we include 
no quark-exchange diagram and hence no Pauli-blocking 
effect among quarks is activated. 

At a glance, one can find that the interaction in this 
channel is always attractive, and that the strength of this 
attraction depends on quark masses. The i?-depcndcncc 
is monotonous at all the k's. In the large R region 
(i? > 4), the potential has smaller m q dependence, while 
it has strong m q dependence in the short-distance region 
(R < 4). The potential V(R) rapidly reduces with de- 
creasing k (increasing m q ), and finally the potentials at 
k=0.1350 and 0.1400 coincide with each other. The po- 
tential exhibits saturation at heavy quark-mass region, 
which implies a long-range quark-mass independent at- 
traction. Although quark masses are also responsible for 



D. Potentials : (i, j, k, I) = (1, 2, 1, 2) case 

The interhadron potentials with the flavor combina- 
tion (i,j,k 7 l) = (1,2,1,2) can be seen in Fig. [5] The 
difference from the previous (i, j, k, I) = (1, 2, 3, 4) case is 
that quark-exchange diagrams are now included, which 
gives rise to Pauli-blocking effect among quarks. 

We readily find that strong repulsions in the short- 
distance region appear in this case. Compared with 
the (i,j,k,l) = (1,2,3,4) case, the quark-mass depen- 
dence is not monotonous: They are smooth functions of 
R at large k's, while the potential at k=0.1350 has a 
pocket at the intermediate distance. As the quark mass 
decreases, the short-range repulsive part rapidly grows 
up and the intermediate attractive pocket disappears. 
Such qualitative change in shape suggests that the po- 
tentials in (i,j,k,l) = (1,2,1,2) case consist of two or 
more parts; attractive part and (probably) repulsive part 
whose strengths and quark-mass dependences are differ- 
ent from each other. 

E. Potentials : Quark-exchange part 

In the following, we define a "quark-exchange part" 
V cx (R,m q ) in V 12 s 2 (R,m q ) as 

V ex (R,m q ) = Vi2,w(R,in q ) - Vi2,3±(R, m q ), (5) 

which is nothing but the difference between ^12,12 (R, ™ 9 ) 
and Vi2,34(-R, rn q ). Taking into account that quark- 
exchange diagrams are included only in the case 
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FIG. 6: Potentials Vi2,i2(-R) computed with the flavor com- 
bination, (i,j,k,l) — (1,2,1,2), are plotted as functions of 
relative distance R. 

°f Vi2,i2(-Rj m q) and that the diagrams needed for 
V12, 34(R, nT>q) are identical with the direct diagrams in 
^I2,i2(-Rj m q)i it gives a measure of "quark-exchange ef- 
fect" or Pauli-blocking effect among quarks in the poten- 
tial Vi2,i2(-R, trig). In practice, it is equivalent to assum- 
ing 

V r i2,34(-R, m q ) = Vdir(-R, m q ), (6) 

V12, 12 (R,m q ) = Vdiv{R,rn g ) + V ex (R,m q ), (7) 

for Vi2 t 34(R,m q ) and Vi 2t i 2 (R, m q ), with a "direct 
part" Vdi r (i?, m q ) measured only with direct diagrams 
and a "quark-exchange part" V CX (R, m q ) induced by 
adding exchange diagrams or equivalently by introduc- 
ing Pauli blocking among quarks. Such decomposition 
of a diquark-diquark potential is conceptually similar 
to the usual decomposition procedure in RGM (Res- 
onating Group Method) calculations in quark cluster 
models [ll|, [H, |3l|. In case multiply repeated quark- 
exchange processes can be neglected, like the Born ap- 
proximation, these two components, and V CX1 result 
in so-called direct and exchange potentials, respectively. 
Thus extracted potential V ex (R,m q ) = Vi2,i2(R,m q ) — 
Vi2,34(R, TUq) shows a monotonous behavior as is seen in 
Fig. and therefore the decomposition seems beneficial 
for our purpose to clarify strengths or ranges of poten- 
tials. At the same time, one finds that the short-range 
repulsion arises only from V ex (R,m q ). Pauli-blocking ef- 
fect among quarks is essential for the short-range repul- 
sion. 



IV. DISCUSSIONS 
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FIG. 7: "Quark-exchange parts" of potentials, which are de- 
fined as V ex (R) = Vi2,i2(-R) — Vi2,34(-R), are plotted as func- 
tions of relative distance R. 



expressed as 

Vi2,u(R, m q) = Vdii(i2, m q ), (8) 
V 12 ,i2 {R, m q ) = V d i I (R,m q ) + V ex {R,m q ). (9) 

Vdi T (R, "m q ) and V ex (R, m q ) can be regarded as the origins 
of attractive and repulsive interactions, respectively. 

The key property of the independent potentials is 
twofold; interaction ranges and strengths. Their m q - 
dependences are also helpful to clarify the origins of the 
potentials. In phenomcnological models, hadronic inter- 
actions are usually incorporated as one-boson-exchangc- 
potcntials (OBEP), which are Yukawa-type and whose 
interaction ranges are directly connected to (exchanged) 
meson masses. Then, such OBE potentials should be 
sensitive to meson masses, i.e. quark masses. On the 
other hand, the short-range repulsive core is often consid- 
ered to arise from a color-magnetic interaction induced 
by one-gluon-exchangcs (OGE), whose strength is pro- 
portional to the inverse square of (probably constituent) 
quark masses, ~ to q 2 - Such ?7iq-dependent interactions 
could be accessed by monitoring the strengths of the po- 
tentials. 

In this section, we evaluate strengths and ranges 
of these decomposed potentials, Vdi r (-R, rn q ) and 
V ex (R,m q ), together with their m g -dependences. Since 
we do not have priori function forms, especially for the 
short-range repulsive part, we try the following five trial 
functions, Fi(x) (i = 1, 2, 3, 4, 5), which have two param- 
eters, strength A and range B. 



Now that we have extracted the ingredients of each 
potential. Namely, Vi2,i2(R, m q ) and Vi2.34(R,m q ) are 
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The fit results can be all found in Table. [TT1 in which 
X 2 /7Vdf, strengths A and ranges B for Vdi r (-R, m q ) and 
V^ x (i?,m 9 ) at each n are listed. None of five functions 
produces x 2 /Wdf — 1 for all the data, which would be 
due to the lack of knowledge about detailed function 
forms, some systematic errors in the lattice data, and 
possible multi-range components, which will be discussed 
later. We adopt F3 as the fit function, since it yields rel- 
atively smaller x 2 /./Vdf for any combination of potential 
type and quark mass. For reference, we plot (best-fit) 
i*3(ii)'s and the lattice data in Fig. [8] The lattice data 
are rather well mimicked by the function. 



A. strengths and ranges : Vdh 




relative distance [lattice unit] 



In Table |TT] and Fig. |H1 we show the fitted parame- 
ters, strength A^h and range Bdir of the attractive part 
Vdir(-R, TTLq), as functions of half of axialvector-diquark 
mass. (Shown as "DIR" in Fig. O) Interestingly enough, 
both of the strength and the range of Vdi T (R, m q ) ex- 
hibit flattening in the heavy quark-mass region, which 
indicates that Vdi r (-R, m q ) contains a universal attrac- 
tive potential V^ t (R). "Universal" here means that nei- 
ther the strength nor the range of V a u tt {R) depends on 
quark mass and that V att {R) always appears in any fla- 
vor channels. In fact, an attractive "dip" at the inter- 
mediate range in Vi2,i2(-R) Tn q ) seems to originate from 
this universal attractive potential V a u tt (R). The combina- 
tion of TOg-depcndcnt short-range repulsive force and m q - 
independent long-range attractive force provides a dip at 
the intermediate range. 

In much heavier quark-mass region, the properties of 
quark wavefunctions would undergo gradual change, and 
V a u tt {R) could be quark-mass dependent. However, in 
the quark- mass range we consider here, V att {R) can be 
treated as a fixed function. 

Then, the next question is what the reminder is. We 
define the quark- mass dependent part V a ® t (R, m q ) as 

V a Z(R, m q ) = V dil (R, m q ) - V&(R). (15) 

We also fit V a ® t (R,m q ) with F3, and the fitted parame- 
ters are shown in Table [TT] and Fig. |H1 Here, we simply 



FIG. 8: Upper. Attractive parts of potentials, Vdir(-R> Wj) = 
Vi2,34(R,m q ), together with the best-fit curves Fs(R) are 
plotted as functions of relative distance R. Lower. "Quark- 
exchange parts" of potentials, V C x(R,m q ) — Vi2,i2(R,m q ) — 
Vi2,34{R,rn q ), together with the best-fit curves are plotted as 
functions of relative distance R. 

adopt Vdi r (R,m q ) at k — 0.1350 as the universal part 
V^ t (R), since we observe almost no quark-mass depen- 
dence already at this k. The strength and the range 
of V£ t {R,m q ) are shown as "ATT(D)" . The extracted 
range B att is no longer dependent on quark mass. That 
is, Vdir(R,m q ) is approximately described by two inde- 
pendent parts; 

V diT (R,m q ) ~ VV t (R) + V a °(R,m q ) (16) 
= AljV t {R)+A° t {m q )f° t {R). (17) 

First part represents a m g -indepcndent weak long- 
range force, and the second one does a short-range re- 
pulsive interaction that has an m 9 -dependent strength 
and an m g -insensitive interaction range. Possible ori- 
gins of the universal attraction would be transition pro- 
cesses to other two hadronic (intermediate) states, which 
is schematically illustrated in the left panel of Fig. [TU1 or 
ra g -independent gluonic interactions. As the candidates 
for m 9 -indcpcndcnt gluonic interactions, one can con- 
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V dir 


0.1350 


0.1400 


0.1450 


0.1500 


X 2 /N d{ (F 2 ) 
X 2 /N d f(F 3 ) 
X 2 /NMF4) 
X 2 /N df (F 5 ) 


1.462 
54.50 


1.719 
6.970 


12.24 
0.381 


3.163 
9.037 


8.938 
33.78 
85.32 


1.094 
3.725 
15.00 


3.516 
0.094 
3.818 


2.206 
11.33 
30.43 


MFi) 
A(F 2 ) 
A(F 3 ) 
A{F A ) 
A(F 5 ) 


0.046(001) 0.050(002) 0.081(009) 0.345(007) 
0.015(004) 0.017(003) 0.032(001) 0.144(010) 
0.154(010) 0.171(006) 0.294(019) 1.134(024) 
0.073(013) 0.083(008) 0.162(002) 0.655(050) 
0.011(008) 0.015(007) 0.046(007) 0.206(052) 


B(F 1 ) 

B(F 3 ) 
B(F 4 ) 
B(F 5 ) 


4.24(006) 
7.06(150) 
1.93(008) 
3.27(041) 
10.3(060) 


4.33(011) 
6.77(076) 
1.91(005) 
3.11(021) 
8.12(277) 


3.87(024) 
5.36(010) 
1.62(007) 
2.42(002) 
4.63(050) 


2.58(005) 
3.61(020) 
1.16(002) 
1.72(011) 
3.20(058) 




0.1350 


0.1400 


0.1450 


0.1500 


X 2 /N df (F2) 
X 2 /N d i(F 3 ) 
X 2 /N df (F 4 ) 
X 2 /N di (Fs) 


3.869 
0.546 


5.317 
1.144 


6.064 
3.713 


7.543 
4.562 


0.095 
1.871 
8.460 


0.481 
4.896 
16.28 


0.737 
8.627 
27.94 


2.202 
13.19 
33.92 



A(Fi) 
A{F 2 ) 
A(F 3 ) 
A{F A ) 
A(F 5 ) 


0.060(003) 0.099(006) 0.209(012) 0.565(036) 
0.025(001) 0.046(002) 0.094(006) 0.277(016) 
0.201(002) 0.354(008) 0.744(017) 2.188(085) 
0.119(007) 0.223(020) 0.453(046) 1.424(163) 
0.040(010) 0.089(023) 0.168(050) 0.606(152) 


B{F X ) 
B(F 2 ) 
B(F 3 ) 
B(F t ) 
B(F 5 ) 


2.58(010) 2.78(010) 2.99(010) 3.25(007) 
3.58(008) 3.63(009) 3.95(015) 3.97(009) 
1.14(001) 1.17(002) 1.27(002) 1.29(002) 
1.67(008) 1.64(009) 1.80(012) 1.72(008) 
2.98(048) 2.61(041) 2.98(055) 2.60(030) 


''att 


0.1350 0.1400 0.1450 0.1500 


xVNdi(F 3 ) 


N/A N/A 1.338 0.246 


A{F 3 ) 


N/A N/A 0.186(024) 1.041(046) 


B(F 3 ) 


N/A N/A 1.05(009) 1.01(004) 



TABLE II: The best-fit parameters for V ex {R,m q ), 
Vdii(R,m q ), and V® t (R, m q ). A(Fi) and B(Fi) denote the 
strength and the range estimated with the function Fi. 



sider color-Coulomb interaction, medium-range attrac- 
tions due to flavor exchange processes, and so on, since 
no quark-mass dependence is observed. As is found also 
in our analyses, such a weak universal attractive force is 
readily masked by the strong repulsive force in the Vi2,i2 
potential for lighter quark-mass. 



B. strengths and ranges : Kx 

In Table QI] and Fig. [9j we show the htted parameters, 
strength A cx and range B cx of the quark-exchange part 
V cx {R 1 m q ) 1 as functions of half of axialvector-diquark 
mass. (Shown as "EX" in Fig. |H1) Let us first have a 
look at the range parameter B CXl which are shown as 
"EX" in Fig. |H1 One can find that the range B cx in 
V CX (R, m q ) is almost quark- mass independent. It is inter- 
esting since the repulsive cores are sometimes described 
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FIG. 9: Upper: The fitted strengths, A dil , A cx , and A® u , 
of the potentials, Vdir(R, m q ), V cx (R,m q ), and V^ t (R, m q ), 
are plotted as functions of half of axialvector-diquark mass. 
Lower. The fitted interaction ranges, Sdir, Box, and B® tt , of 
the potentials, Vdii(R,m q ), V ex (R,in q ), and V att (R, m q ). The 
parameters for Va\v(R, rn q ), V ex (R,m q ), and V a , it {R, m q ) are 
respectively shown as "DIR", "EX", and "ATT(D)". 



by e.g. omega-meson exchanges in phenomenological 
models. The constituent quark- mass (half of AV-diquark 
mass) variation from 0.53 to 0.22 gives rise to only 10% 
change in the range B exi which is much smaller than ex- 
pected from the meson-mass change. Moreover, at two 
largest k's (two lightest quark masses), the range B ex 
remains unchanged. V cx (R,m q ) can be predominantly 
expressed as 



Ve X (i?, m q ) - A ex (m q )f ex (R). 



(18) 



The quark-mass dependence (approximately) appears 
only in the strength A cx . From these observations, we 
conjecture that the short-range repulsion is not generated 
by meson poles as OBEP, but by some other origins. 

One most possible mechanism for the short-range re- 
pulsion is a color-magnetic (CM) interaction among 
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quarks as suggested by constituent-quark models. A CM 
interaction accompanied by Pauli blocking among quarks 
raises the energy of (closely located) two hadrons, which 
results in strong repulsions. In the first-order perturba- 
tion, the contributions from the CM interaction are pro- 
portional to its strength, whose m 9 -dependence is given 
by m Q 2 . 

The strengths A ex extracted by fits are also shown in 
Table Ulland Fig. El Acx increases in the light-quark-mass 
region and it is qualitatively consistent with the CM in- 
teraction as the origin of repulsion. In the upper panel 
in Fig. HI two fit functions, Cttiq 2 and Ctoq 3 are respec- 
tively plotted as a long-dashed and a dashed line. The 
strength A cx is well reproduced by Cmg 3 rather than 

Cttiq 2 , which implies that the quark- mass dependence 
of repulsion seems stronger than that of the strength of 
the CM interaction itself. However, considering that the 
Pauli blocking among quarks is essential for the repulsion 
and that the interaction range of the repulsive force has 
a very weak quark-mass dependence, the present results 
are consistent with the quark-model interpretation that 
short-range repulsion between hadrons arises from CM 
interaction among quarks. 



V. SOME ADDITIONAL TRIALS 

A. Boundary condition dependences 

In the previous section, we claimed that one possi- 
ble contribution for the universal attraction would be 
transition processes to other two hadronic (intermediate) 
states, which is schematically illustrated in the left panel 
of Fig. [TO] In case that two hadrons are closely located, it 
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FIG. 10: Schematic figure for quark-exchange interactions. 
Two hadrons, whose flavors are initially and finally (1, 2) and 
(3, 4), interact with each other through the intermediate state, 
where their flavors are (1,3) and (2,4), for example. If we 
impose the antiperiodic boundary condition onto 3- and 4- 
quarks, (1, 2)- and (3, 4)-diquarks obey the periodic boundary 
condition, whereas the intermediate (1, 3)- and (2, 4)-diquarks 
do the antiperiodic boundary condition. As a result, the en- 
ergies of the intermediate states are raised up, if they are 
spatially spread and being affected by finite-size effects. 

may not be meaningful to identify their flavors, because 
in such short-distance regions, two hadrons are largely 



overlapping with each other and they may differ from 
two isolated hadrons [IH HU ■ 

In order to get sure that we have no finite- volume arti- 
fact for the universal attraction, we repeat the same anal- 
yses with antiperiodic spatial boundary conditions [HI, 
|34| |. The flavors arc again (i,j,k,l) = (1,2,3,4), and 
anti-periodic boundary conditions are imposed for 3,4- 
flavor quark fields. In this case, (l,2)-(3,4) two-diquark 
state is never affected, since (1,2)- and (3, 4)-diquarks 
obey a periodic boundary condition. If there exists finite- 
volume artifact in intermediate transition states, e.g. 
(1, 3)-(2, 4) hadronic state, the results would be modified 
in this new analysis. We eventually observed no modifi- 
cation in all the cases, and no serious finite volume effect 
is confirmed. 



B. Removal of high- momentum gluons 

Hadronic interactions may be classified into two cate- 
gories. One is highly nonperturbativc phenomena, such 
as meson interactions, and the other is perturbative con- 
tributions represented as CM interactions among quarks. 
In one sense, these phenomena differs in energy scale 
they belong to. Recently, a method that restricts en- 
ergy scale of gluons via the Fourier transformation was 
proposed and tested [3|1 HI] . The authors reported that 
removal of high-energy contributions results in vanishing 
short-range Coulomb interaction in heavy-quark poten- 
tials, keeping the linear confinement part still unchanged. 
Taking into account that the Coulomb potential in short- 
distance region comes from OGE processes, removal of 
high-energy gluons is expected to cut such OGE pro- 
cesses. Though it should be clarified whether such re- 
moval also cuts other OGE interactions, e.g. color- 
magnetic part, this analyses might give us a hint to 
short-range hadronic interactions. Especially, the ori- 
gin of short-range repulsive core coming from the color- 
magnetic OGE interaction would be accessed to some 
extent in this manner. 

We perform 3-dim Fourier transformation [37| for 
link variables (x, t) in the Landau gauge and extract 
Ufi(p,t) in momentum space: 

u^t) = Y^u^,tyv* (19) 

X 

We leave the low-momentum contribution, 



u low (p,t) = I u ^ l) (|p| - A) 

M ' I (otherwise) 



(20) 



and reconstruct new SU(2) link variables in the coordi- 
nate space, [7^ ow (x, t), so that the distance, 

Tr (£# w (x,i) - C7^(x,t))(t/^(x,t) - ^ ow (x,t))t ! 

(21) 
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is minimized. Here, 

U]T(x, t) = - £ ^ W (P' *) e " P ' X ( 22 ) 
p 

is a truncated link variable reconverted into the coordi- 
nate representation, which no longer belongs to SU(2). 
The infrared cut A is set to A = 5 in lattice unit. 




distance r [lattice unit] 



list the masses of scalar (S) and axialvector(AV) diquarks. 
The masses measured with low-momentum link variables 
almost show no mass splitting, which implies that the 
color-magnetic interaction that gives rise to S-AV mass 
splitting is now largely suppressed. The fact that S- and 
AV-diquark masses degenerate in this manner again im- 
plies the validity of nonrclativistic-quark-modcl descrip- 
tion of the diquarks. The NG-boson nature in S-diquark 
hardly appears in the present quark-mass region. 
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FIG. 11: Heavy quark-antiquark potential as a function of 
a separation r. The upper data shown as circles denote the 
lattice QCD data obtained in full SU(2) calculation, and the 
lower data shown as squares denote those obtained without 
high-momentum gluons. The infrared cut A is set to A = 5 in 
lattice unit. The original QQ potential is fitted as Vqq(t) = 
-0.2397(53) i + 0.0518(10)r + 0.5405(50). 

We show the intcrquark potential in Fig. 1111 The up- 
per data in Fig. Qj] shown as circles denote the lattice 
QCD data Vqq(r) obtained in full SU(2) calculation, and 



the lower data shown as squares denote those V, 



low / 

QQ 



obtained without high-momentum gluons. The original 
QQ-potential data Vqq(r) is fitted with the Cornell type 
potential as, V QQ {r) = -0.2397(53)^ + 0.0518(10)r + 
0.5405(50). On the other hand, a fit of V%%(r) using 
the Cornell type potential seems not applicable due to 
the oscillation behavior in the short-range regions. We 
simply draw a dotted line, V(r) = 0.0518r, in the figure. 
One can find that the short-distance Coulomb interaction 
disappears and the linear confinement part remains al- 
most unchanged. At the same time, the constant part in 
the potential, which represents the self energy of a static 
fundamental charge under the lattice regularization, also 
disappears in this treatment. 

An interesting observation can be found in the masses 
of scalar and axial vector diquarks. We adopt two hopping 
parameters K = 0.1150 and 0.1250. We note here that we 
have adjusted k's to cover approximately the same quark- 
mass region so that the AV-diquark masses before and af- 
ter high-momcntum-gluon cut are similar. In Table. HI we 



FIG. 12: Potentials V12, 34(A) for the flavor combination 
(i, j, k, I) — (1, 2, 3, 4) computed without high-momentum glu- 
ons are plotted as functions of relative distance R. The in- 
frared cut A is set to A = 5 in lattice unit. 

We show the interhadron potentials measured with 
the low-momentum link variables in Figs. I12[ 1131 and 
1141 Naive expectation is that potentials originating from 
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FIG. 13: Potentials Vi2,i2(R) for the flavor combination 
k, I) = (1, 2, 1, 2) computed without high-momentum glu- 
ons are plotted as functions of relative distance R. The in- 
frared cut A is set to A = 5 in lattice unit. 
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taincd without high-momentum gluons. Open and filled 
symbols are compared at the same Am. While the open 
symbols show Am5 ~ m Q 3 behavior, the filled symbols 
deviate from the expected curve and always lie above it, 
which indicates that the repulsive interaction does not 
always decrease like the S-AV mass splitting, and that 
some extra enhancement remains even after the high- 
momcntum-gluon cut. Wc then expect that the repulsive 
force still contains small contributions other than CM in- 
teraction among quarks. 

Restricting gluon's energy scale surely changes poten- 
tials and may be useful to clarify detailed structures. 
However, there seem to remain several issues to be clari- 
fied before quantitative conclusions are reached. 



C. Different choice in operators 



FIG. 14: "Quark-exchange parts" of potentials, which are 
defined as V EX (R) = Vi2,i2{R) — Vi2,34(-R), computed without 
high-momentum gluons are plotted as functions of relative 
distance R. The infrared cut A is set to A = 5 in lattice unit. 

short-range OGE processes are suppressed, while those 
from other origins remain unchanged. Actually, all the 
potentials get smaller after the high-momentum gluons 
are removed, which can be seen in Figs. [I2j HSl and [14] 
The i?-dcpcndence of potentials are also much milder 
now, and the divergent behavior around the origin dis- 
appears. The reason why all are now decreased could 
be that many of the hadronic processes, e.g. flavor ex- 
changes, meson exchanges, and so on, need color recom- 
binations in hadrons. If such recombinations are prohib- 
ited, some of hadronic processes depending on them will 
be also suppressed. 




-0.2 1 1 1 1 1 1 1 1 1 

0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 

Am [lattice unit] 

FIG. 15: The values of V c * at R=l and 2 are plotted as 
functions of S-AV mass splitting Am. 

To clarify the origin of repulsion in a clearer way, we 
plot the values of V oyi at R=l and 2 as functions of S- 
AV mass splitting Am in Fig. \T5\ The open symbols are 
the original data taken from fully calculated potential 
V cx , and the filled symbols denote the values of V cx ob- 



Concern about possible operator dependences of po- 
tentials from BS amplitudes has been often raised [38j j . 
BS amplitudes themselves are operator-dependent quan- 
tities, and inevitably equivalent potentials extracted from 
them are. Though the strengths and the ranges of the 
potentials would be operator dependent, some essential 
features in hadronic interactions will be clarified by inves- 
tigating operator dependences. Among possible hadronic 
operators, we hereby consider spatially extended smeared 
operators, because hadron-sizc effect is considered to be 
one of the most important features: Hadronic states have 
finite sizes, and hence finite-size (smeared) hadronic op- 
erators would be more adequate to describe hadrons. 
We note that the reduction formula needed to relate 
BS amplitudes to scattering observables has been proved 
only for point-type operators [39l - l42| . In this subsection, 
we repeat our analyses employing smeared operators for 
sinks in order to derive essentials in hadronic interactions. 

First, we mention smearing methods. Usually, quark 
operators are individually smeared, and hadronic opera- 
tors are constructed from them. In the case of gauge- 
invariant Gaussian smearing, each quark operator is 
smeared with a gauge-covariant lattice derivative oper- 
ator 

K^{U) XV = U^(x)6 x+ p, ty + U^x - (i)5 x -^y (23) 

as 

g 6 (x;r) = H^l + aJ2 K i(UUyJS y , r . (24) 

Parameters a and are chosen so that quark fields are 
distributed around r with a radius b. We hereby con- 
sider diquark states which consist of two quarks with 
identical masses. A smeared hadronic (diquark) oper- 
ator i?b(xi,X2;x) whose center position is x is then de- 
fined with gi j (,(xi;x) and g2,6(x2;x) as Hb(xi, x 2 ; x) = 
<2i.b(xi; x)<7 2 ,f)(x2; x). However, in reality, iJ b (xi,x 2 ;x) 
does not describe a hadronic state located at x, but it 
contains unwanted contributions because a "position" of 
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a hadron, (xi + X2)/2, docs not always coincide with 
x. For example, i^fxii x 2 ; x) has nonzero entry at 
xi = X2 = x' 7^ x, which is nothing but a point-type 
hadronic operator located at x' 7^ x. While such con- 
tributions are not harmful in hadron-mass measurement 
because hadronic states arc projected onto a momentum 
cigenstate summing up hadronic correlators, it causes se- 
rious difference in the measurement of BS amplitudes. 
-ff&(xi,x 2 ;x) should be zero when (xi + x 2 )/2 7^ x. In 
order to satisfy the condition (xx + x 2 )/2 = x, we im- 
prove smearing function with the following prescription. 

With new link variables, U^ 1 = exp(i0 M )L^, and a 
phase < 8^ < 2tt, we define 



rfrWxv (25) 
= Uy(x)S x+ ^ y + U^(x-fi)S x .^ y (26) 

= exp(-i6n)Ufa)6 x+ ji ty + exp(i0 M )J7 M (a; - (t)6 x 
Smeared quark operators are similarly defined as 

AT / 3 \ 



^(x;r) = J] 1 + «E* 



(28) 



! = 1 



Here, 6 = (61,82,63) can be chosen freely. As a result, 
<^(x; r) are different from gt,(x; r) in overall phases; 



q e b (x; r) = exp(id ■ (x - r))^(x; r) 



For simplicity, we hereby consider a hadronic operator 
located at the origin, x — 0. Then, a hadronic (diquark) 

operator H b l 2 (xi, X2; 0) is now 

^f A ( Xl ,x 2 ;0) (30) 
= g f 6 ( Xl ;0) g |f 6 (x 2 ;0) (31) 
= cxp(i^i ■ xi + id 2 ■ x 2 )qi^(xi; 0)(72,b(x 2 ; 0). (32) 

Taking 8 X = 6 2 = 8, 

fff(x l5 x 2 ;0) (33) 
= exp(i6 ■ (xi + x 2 ))qi^(xi; 0)q 2 .b(x 2 ; 0). (34) 

If we average 7J^(xi,x 2 ;0) as ^^_ff^(xi,x 2 ;0) over 8 

so that exp(i6* • (xi + x 2 )) = <5 X i+x 2 ,o, the "position" 
of a hadronic operator can be unambiguously defined, 
(xi + x 2 = is always satisfied.) The simplest choice 
is averaging it over a random set {8}. This scheme is a 
generalization of projection processes, and could be em- 
ployed for other purposes. 

We note here that nothing other is changed in this 
process. A hadronic operator is now "projected" onto 
(xi + x 2 = 0)-statc so that unwanted contributions are 
absent. In actual calculations, _ff^(xi,x 2 ;0) are aver- 
aged over the set {6} = {jf-(l, m, ri)}(l < l,m,n < N$). 
When Ng is smaller than the lattice spatial extent, the 



delta function (5 X i+x 2 ,o is not fully reproduced and the 
projection is incomplete. To eliminate such (small) con- 
taminations, we further introduce overall random phases. 

Next, we measure hadron "sizes" monitoring overlap 
coefficients. The sink operators relevant to BS ampli- 
tudes are now smeared in a gauge invariant way. The 
root-mean-square radius 6 smr of the smeared operators 
are determined so that the smeared operators have max- 
imal overlaps with the ground states (scalar diquarks). 
In Fig. [THl we plot the squared overlaps, |cq| 2 / l c i| 2 - 




12 3 

RMSR of smeared operators [lattice unit] 



(29) FIG. 16: jc | 2 / J2i M 2 are plotted as a function of the RMSR 



of a smeared operator 6 smr . Here, d 
diquark operator and i-th state; a = 



denotes the overlap of a 
(vacl-Dj^rlith state). 



Here, Cj denotes the overlap of a diquark operator and 
i-th state, which is defined as 



^ = (vacl-D^rKth state). 



(35) 



While the squared overlaps are slightly dependent on 
quark masses, they are saturated at 6 smr ~ 3. We then 
set 6 smr = 3 for all the k's. In Figs. [T7]fT91 we show 
Vi2,34,Vi2,i2, and V cx obtained with "properly projected" 
smeared operators. 

V12.34 obtained with smeared operators are qualita- 
tively similar to those with point operators. It is always 
attractive and gets stronger at smaller quark-mass re- 
gion. Especially, the saturation at heavy tuq region can 
be again seen in Fig. [T7] The interaction range of the 
universal attractive potential V^ t (R) is now reduced by 
about 20%, and the strength is stronger, compared with 
the universal attraction observed in the point-operator 
case. On the other hand, repulsive cores in V\i.\i at 
short-distance region are generally weaker in this case, 
which is displayed in Fig. 1181 Apparent repulsive cores 
disappear at small k's, whereas cores still persist and are 
observed at larger k's. V\2,i2 seems more sensitive to 
operator choice than ^12,34. Interestingly, the difference 
V ex (R) = Vi2,i2(-ft) — ^12,34 (-R) still hold similar prop- 
erties to that in point-operator case. It grows as tuq 
decreases and is always repulsive, which confirms that 
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the Pauli-blocking effect is responsible for repulsion also 
in this case. The interaction range of V cx remains almost 
unchanged as compared to that extracted with point op- 
erators. 

Quantitative evaluation of interhadron potentials 
seems difficult at present, because of its operator de- 
pendences. However, several basic properties remain 
qualitatively unchanged with/without operator smear- 
ing: (1) The existence of universal attraction between 
two hadrons. (2) The difference V ex (R) = Vi2.i2{R) — 
Vi2,3i(R) still shows repulsive contribution at any k's. 

(3) The quark-mass dependences of the interaction ranges 
of direct (Vi2,34) and exchange (V ex ) parts of potentials 
are much weaker than expected from meson exchanges. 

(4) The quark-mass dependence of the strength of V cx is 
consistent with or stronger than m^ 2 . 
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FIG. 18: Potentials Vi2,i2(-R) computed with "projected" 
smeared diquark operators ^2gH®, with the flavor combina- 




FIG. 17: Potentials V12, 34 (R) computed with "projected" -0.2 1 1 1 1 1 1 1 1 1 1 

smeared diquark operators J2gHb> witn the flavor combina- relative distance [lattice unit] 

tion, k, I) = (1, 2, 3, 4), are plotted as functions of relative 
distance R. 



VI. SUMMARY AND OUTLOOKS 

We have evaluated inter-hadron interactions in two- 
color lattice QCD based on the independent function 
V(R). The function V(R), which has been expressed as 
"potential" throughout this paper, has been extracted 
from the Bethc-Salpctcr wave functions on the lattice as- 
suming a nonrclativistic Schrddinger equation. The simu- 
lations have been performed in quenched two-color QCD 
with the plaquette gauge action at f3 = 2.45 and the 
Wilson quark action. Evaluating different flavor combi- 
nations as well as several quark masses, we have extracted 
and investigated each ingredient in hadronic interactions. 
We have considered two-diquark scatterings, whose fla- 
vors are and (k,l). When all the quarks have dif- 
ferent flavors, (i, j, k, I) = (1, 2, 3, 4), the interhadron po- 



FIG. 19: "Quark-exchange parts" of potentials, which are 
defined as Kx(fi) = Vi2,i2(R) — Vi2,34{R), are plotted as 
functions of relative distance R. In this case, Vi2,i2(R) and 
Vi2,34(R) are those measured with "projected" smeared di- 
quark operators ^2gH®. 



tential is quark-mass dependent and always attractive, 

Vl2,34 (R, m q ) = Vdir(-R,TOq). (36) 

On the other hand, in (i,j,k,l) = (1,2,1,2) chan- 
nel, where exchange diagrams, i.e. Pauli-blocking ef- 
fect among quarks, are included, Vi2,i2(-R, m q ) shows 
repulsion in short-distance region, and this repulsion 
arises only from the exchange diagrams. We have de- 
fined the "quark-exchange part" V CX (R, m q ) in the po- 
tential, which is induced by adding quark-exchange di- 
agrams, or equivalcntly, by introducing Pauli blocking 
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tributions can be further decomposed into simpler parts; 
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FIG. 20: Upper. The fitted strengths, A dil , A cx , and A£ u , 
of the potentials, Vdi r {R, m g ), V ex (R,m q ), and V£ t (R,m g ), 
are plotted as functions of half of axialvector-diquark mass. 
Lower. The fitted interaction ranges, Bdir, B cx , and B® tt , of 
the potentials, Vdi?(R, m q ), V ex (R,m q ), and V aU (R, m q ). The 
parameters for Vdir(R, m q ), V ex (R,m q ), and V^i t (R,m q ) are 
respectively shown as "DIR", "EX", and "ATT(D)". 
All the parameters are obtained by fitting potentials measured 
with "projected" smeared diquark operators ^2gH® • 



among some of quarks. Namely, when some of the quarks 
have identical flavors, (i,j,k,l) = (1,2,1,2), the quark- 
exchange part V cx (R,m q ) is defined as V ox (R,m q ) = 
Vi2,i2(R, m q ) — Vi2,3i(R, Trig), which means the inter- 
hadron potential is written as 



Vi2,i2(R,m q ) = Vdir(R,rn q ) + V cx (R,m q ). 



(37) 



It indicates that the short-range repulsive core in 
Vii,i2 (R, TUq) arises only from the "quark- exchange part" 
V cx (R,m q ). Pauli blocking among quarks seems essen- 
tially needed to generate repulsive force between two 
hadrons. 

We have found that these R- and m 9 -dcpendent con- 



V ex (R,m q ) ~ A eyL (m q )f ex (R), (38) 
V dlI (R,m q ) ~ VV t (R) + V£(R,m q ) (39) 
= AVjg t (R)+Ag t (m g )f° t (R). (40) 

Here, f's are quark-mass independent functions of R. In- 
terestingly, Vdir(i?, m q ) can be further decomposed into 
Vf/ tt (R) and Vj^ t (R, m q ), a universal attractive part and 
an m g -dcpcndent attractive part. 

One prominent observation is that all the JJ-functions 
fex(R), f^ tt (R), and f^ t (R) arc quark-mass indepen- 
dent, which implies that the quark-mass dependences of 
V ex (R,m q ) and Vj^ t (R,m q ) do not appear in interaction 
ranges but only in the overall strengths. Though we first 
adopted the specific function F^(x) for potential- form 
analyses, we finally found that simple rescaling of the po- 
tential V ex (R,m a ) (the lattice data themselves) well ex- 
plains all the quark-mass dependences. That is, we even- 
tually encounter no scheme dependence in these descrip- 
tions. These observations apparently do not go together 
with the meson-exchange picture among hadrons, since 
TO 9 -dependences of interaction ranges seem very small. 
At least in the present quark-mass range, mesonic con- 
tributions seem small and subdominant, and hence we 
observed interactions of the other origins. 

For the repulsive part, we have found that the strength 
grows as ~ m q 3 > which is similar to but slightly stronger 
than the color-magnetic interaction itself by onc-gluon- 
exchange (OGE) processes. Considering that the Pauli 
blocking among quarks is essential for the repulsion and 
the interaction range of the repulsive force has a very 
weak quark-mass dependence, it is likely to originate 
from a color-magnetic interaction among quarks. 

Actually, the function V(R) is energy and operator 
dependent. In order to clarify operator dependence, 
we have constructed "projected" smeared operators and 
compared the results with and without operator smear- 
ing. As a result, we observed quantitative changes in 
potential shapes. Quantitative evaluation of interhadron 
potentials seems difficult at present, because of its op- 
erator dependences. On the other hand, several basic 
properties remain qualitatively unchanged with/without 
operator smearing: (1) The existence of universal attrac- 
tion between two hadrons. (2) The difference V CX (R) = 
Vi2,i2(-R) — Vi2.3i{R) still shows repulsive contribution at 
any k 's. (3) The quark-mass dependences of the inter- 
action ranges of direct and exchange parts of potentials 
(Vdir,V ex _) are much weaker than expected from meson ex- 
changes. (4) The quark-mass dependence of the strength 
of V cx is consistent with or stronger than tuq 2 . To deter- 
mine the precise forms of interhadron potentials, it would 
be necessary to perform more sophisticated analyses tak- 
ing into account energy dependence in a less operator- 
dependent way. 

Since quark-mass dependence in interaction ranges 
does not appear, the origin of the interactions we ob- 
served so far would be all gluonic interactions and/or 
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flavor-exchange processes, e.g. color-magnetic or color- 
Coulomb interactions, and so on, rather than mesonic 
contributions. As was found in our analyses, attractive 
forces are readily masked by the strong repulsive force 
appearing in the V12.12 potential in the light quark-mass 
region. If the universal attraction in hadronic interaction 
appears also in SU(3) QCD, they might be observed in 
a channel with no or less quark-exchange contribution 
between hadrons, such as Nip scattering state. In the 
lighter quark-mass region, meson-exchange contributions 
could largely emerge and be predominant, which is left 
for further studies. 
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